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Abstract 

A  previously  developed  line-by-line  spectral 
model  has  been  extended  to  use  the  HITRAN  and 
HITEMP  line  strength  data  bases.  The  new  spec¬ 
tral  model  computes  the  line-by-line  spectral  radi¬ 
ance  and  transmittance  of  homogeneous  or  inho¬ 
mogeneous  gas  paths,  including  the  effects  of 
Doppler  shifts  due  to  directed  gas  velocity  and 
vibrational  nonequilibrium.  The  nonequilibrium 
vibrational  state  populations  can  either  be  speci¬ 
fied  directly  (perhaps  by  a  state-specific  chemistry 
model,  or  by  imposing  rotational  and  vibrational 
temperatures  for  each  molecular  constituent).  Cal¬ 
culations  are  compared  to  benchmark  laboratory 
data  for  hot  and  hot-through-cold  gas  paths  of  C02 
at  high  and  low  spectral  resolution.  The  model 
showed  good  agreement  with  the  data,  particularly 
near  the  line  centers.  Some  discrepancies  are 
noted  in  the  line  wings,  perhaps  indicating  a  more 
detailed  model  of  the  line  wings  is  needed.  Sam¬ 
ple  calculations  demonstrating  the  effects  of 
directed  Doppler  shifts  and  vibrational  excitation 
behind  a  normal  shock  are  presented. 

Introduction 

Modeling  molecular  vibrational  relaxation  has 
increased  in  sophistication  in  recent  years.  Tech¬ 
niques  have  evolved  from  simple  phenomenologi¬ 
cal  models  to  multi-temperature  techniques,  to 
state-to-state  chemistry  models.1'3  These  state-to- 
state  chemistry  models  attempt  to  describe  the 
relaxation  process  by  treating  each  vibrational 
state  as  a  separate  species  and  writing  the  appro¬ 
priate  rate  equations  for  the  various  energy 


exchange  processes  (e.g.,  V-T  or  V-V.)  In  this 
manner,  a  distribution  of  states  that  is  initially 
described  by  a  Boltzmann  distribution  will  depart 
from  thermodynamic  equilibrium  when  the  collision 
frequency  is  sufficiently  low.  This  is  in  contrast  to 
multitemperature  models,  wherein  the  vibrational 
state  distribution  is  assumed  to  be  described  by  a 
Boltzmann  distribution  at  one  or  more  tempera¬ 
tures,  while  the  rotational  states  are  in  a  Boltzmann 
distribution  at  another  temperature.  State-to-state 
models  will  eventually  serve  to  validate  or  disprove 
the  multitemperature  approach. 

Computing  the  radiation  emitted  from  these 
non-Boltzmann  distributions  requires  a  radiative 
transfer  model  that  uses  the  vibrational/rotational 
state  populations  as  input.  Such  a  model  was 
developed  previously  in  Ref.  1 .  In  the  current  work, 
this  line-by-line  spectral  model1  has  been 
extended  to  use  the  HITRAN4  and  HITEMP5 
molecular  line  parameter  data  bases.  The  new 
spectral  model  computes  the  line-by-line  spectral 
radiance  and  transmittance  of  homogeneous  or 
inhomogeneous  gas  paths,  including  the  effects  of 
vibrational  nonequilibrium  and  Doppler  shifts  due 
to  directed  gas  velocity.  The  nonequilibrium  vibra¬ 
tional  state  populations  can  either  be  specified 
directly  through  the  upper  and  lower  state  popula¬ 
tions  or  by  specification  of  rotational  and  vibrational 
temperatures  for  each  molecular  constituent. 

Model  Formulation 

The  HITRAN  and  HITEMP  line  strength  data 
files  include  information  on  individual  vibrational/ 
rotational  transitions  for  32  diatomic  and  poly- 
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atomic  molecules.  This  information  includes  spe¬ 
cies  and  isotope  identification,  line  center  fre¬ 
quency,  upper  and  lower  quantum  numbers  and 
branch  assignment,  line  strength  (cm^^molecule/ 
cm2)),  and  transition  moments  (R2  (Debye2)).  As 
shown  in  Ref.  1 ,  radiation  from  a  nonuniform  NLTE 
gas  at  a  particular  spectral  position  is  given  by: 


L  -J[kem(x')-kabS(x')]dx' 

fjkem(x)e  0  dx  +  N0 

o 

l  (1) 

-J[kabs(x)-kem(x)]dx 


where 

f  =  2  x  10~7hc2a>3 


N  is  the  emitted  radiance  in  watts/[cm2  sr  cm-1],  N0 
is  the  incident  radiance,  the  optical  frequency  to  is 
in  units  of  wavenumbers  (cm-1),  Planck’s  constant 
(h)  has  units  of  erg-sec,  the  speed  of  light  (c)  has 
units  of  cm/sec,  and  kabs  and  kem  (described 
below)  have  units  of  cm-1. 


The  absorption  coefficient  is  given  by: 

3  —36 

kabs(x)  =  8iM0  XAicooidiNLi(x)Rfgi(co,  x)  (2) 

where  the  sum  runs  over  all  spectral  lines.  Aj  is  the 
isotopic  abundance  of  the  ith  line,  NLi  is  the  popula¬ 
tion  density  (#/cm3)  of  the  lower  nondegenerate 
rotational-vibrational  state,  R2  is  the  normalized 
transition  moment  listed  in  the  HITRAN  data  base 
in  units  of  Debye2,  gi(co)  is  the  normalized  line 
shape  function  in  units  of  cm  with  a  unit  integral, 
and  d|  is  the  degeneracy  of  the  lower  rotational- 
vibrational  state.  The  line  shape  function  (Voigt, 
Lorentz,  etc.)  carries  an  x  dependency  through 
temperature-  and  pressure-dependent  half  widths. 
The  Voigt  line  shape  approximation  of  Humlicek6  is 
used  in  the  current  model. 


The  stimulated  emission  coefficient  is  given  by: 

3  -36 

kem(x)  =  87t31h°-lAico0idiNUi(x)Rfgi(o),  x)  (3) 
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where  Nyi  is  the  population  density  of  the  upper 
nondegenerate  rotational-vibrational  state  and  dj  is 
the  lower  state  degeneracy  factor.  Vibrational/rota¬ 
tional  nonequilibrium  is  accounted  for  through  the 
state  populations  Nyj  and  NLi.  Doppler  shifts  due  to 
directed  velocity  components  are  accounted  for  by 
shifting  the  line  center  frequencies  (cooi). 

Numerical  Method 

Equation  (1)  is  the  solution  of  the  differential 
radiative  transfer  equation 

dN_N(kem-kabs)  =  fkem  (4) 

dx 

If  the  gas  path  is  discretized  into  n  homoge¬ 
neous  zones,  Eq.  (4)  may  be  solved  numerically  by 


where  8j  is  the  length  of  the  ith  zone.  The  summa¬ 
tion  is  taken  over  the  number  of  zones. 

Benchmark  Computations 

Carbon  dioxide  is  the  only  infrared  active  spe¬ 
cies  considered  in  this  paper.  For  all  the  calcula¬ 
tions  presented  below,  only  the  two  most  common 
isotopes  of  C02  are  considered:  16012C160 
(denoted  626)  and  16013C160  (denoted  636).  The 
terrestrial  abundances  of  these  two  isotopes  are 
0.984260  and  0.01 106,  respectively.  There  are  127 
C02  vibrational  states  included  in  the  HITEMP  data 
base,  with  a  maximum  rotational  quantum  number 
(J)  of  185.  Vibrational  state  populations  for  these 
127  states  were  generated  from  either  thermody¬ 
namic  equilibrium  considerations,  or  by  a  two-tem¬ 
perature  model  in  which  the  vibrational  and  rota¬ 
tional  temperatures  may  be  assigned  different  val¬ 
ues.  No  state-to-state  chemistry  models  were  used. 
This  capability  will  be  demonstrated  in  a  later  paper. 

Comparisons  with  ERIM  Data 

Lindquist,  et  al.  of  the  Environmental  Research 
Institute  of  Michigan  (ERIM)  measured  the  spectral 
radiance  of  a  hot  gas  cell,  both  unattenuated  and 
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attenuated  by  a  cold  gas  cell.7  The  conditions  of 
the  hot  cell  were  1200  K  temperature,  0.10-atm 
total  pressure,  0.01 -atm  C02  partial  pressure,  and 
a  0.60-m  path  length.  The  conditions  of  the  cold 
cell  were  297  K  temperature,  0.13-atm  total  pres¬ 
sure,  0.1 15-atm  C02  partial  pressure,  and  a  100-m 
path  length. 

Comparisons  of  the  ERIM  hot  cell  4.3-p.m  band 
data  with  the  line-by-line  model  calculations  are 
shown  in  Fig.  1 .  The  spectral  resolution  of  the  ERIM 
data  is  approximately  4  cm-1.  The  highly  resolved 
line-by-line  calculation  results  have  been  con¬ 
volved  with  a  triangular  slit  function  to  attempt  to 
match  the  spectral  resolution  of  the  data.  The  gen¬ 
eral  level  of  the  data  and  the  calculations  compares 
favorably,  with  the  most  disagreement  at  lower 
wavenumbers.  The  shape  of  the  computed  band 
differs  significantly  from  the  data  in  the  2275  cm-1 
region.  This  may  be  due  to  either  missing  vibra¬ 
tional/rotational  states  or  transitions  in  the  model. 
The  model  allows  the  influence  of  each  line  to 
extend  to  256  half-widths  (a  typical  half-width  is 
0.07  cm-1).  This  may  be  insufficient  in  some 
regions.  The  disagreement  in  the  low  wavenumber 
wing  of  the  band  is  most  likely  due  to  this  treatment 
of  the  long  tail  of  each  individual  line  or  missing 
vibrational  bands. 


Fig.  1.  Comparison  of  ERIM  hot  cell  data  with 
model  results. 


Comparisons  of  the  ERIM  hot-through-cold  cell 
4.3-pm  band  data  with  the  line-by-line  model  calcu¬ 
lations  are  shown  in  Fig.  2.  The  model  correctly 
computes  the  location  of  the  typical  blue  and  red 
spikes  at  approximately  2400  and  2195  cm-1, 
respectively.  The  blue  spike  radiance  calculations 
are  about  a  factor  of  6  above  the  data.  It  should  be 


noted  that  the  signal-to-noise  ratio  of  the  data  is 
very  poor  in  this  region.  The  disagreement  in  the 
computed  signal  levels  in  the  red  spike  region  is 
probably  due  to  the  treatment  of  the  line  tails  or 
missing  vibrational  bands. 


Fig.  2.  Comparison  of  ERIM  hot-through-cold  cell 
data  with  model  results. 


Comparisons  with  AEDC  Data 

Phillips8  measured  the  high-temperature  trans¬ 
mittance  of  C02  gas  paths  at  a  spectral  resolution 
of  approximately  0.016  cm-1.  These  experiments 
were  conducted  at  the  U.S.  Air  Force  Arnold  Engi¬ 
neering  Development  Center  (AEDC).  The  experi¬ 
mental  conditions  were  1000  K  temperature, 
0.013-atm  total  pressure,  0.001 3-atm  C02  partial 
pressure,  and  a  10.6-cm  path  length.  The  line-by¬ 
line  calculation  results  have  been  convolved  with  a 
triangular  slit  function  to  attempt  to  match  the  spec¬ 
tral  resolution  of  the  data.  A  comparison  of  the 
model  results  and  the  AEDC  data  in  the  2360  cm-1 
spectral  region  is  shown  in  Fig.  3.  The  computed 
transmittance  agrees  very  well  with  the  data  near 
the  line  centers,  but  the  agreement  is  not  as  good 


Fig.  3.  Comparison  of  AEDC  hot  cell  data  with 
model  results. 
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in  the  wings.  This  overprediction  of  transmittance  is 
probably  due  to  the  treatment  of  the  line  tails 
described  above.  The  direct  extension  of  the 
region  of  influence  of  each  spectral  line  is 
extremely  computationally  expensive.  This  is  why 
models  such  as  the  Fast  Atmospheric  Signature 
Code  (FASCODE)9  use  an  approximate  treatment 
of  the  wings.  It  should  also  be  noted  that  the  mea¬ 
surement  of  these  very  high  transmittances  is  diffi¬ 
cult  due  to  source  intensity  drift  with  time  as  well  as 
other  factors,  leading  to  relatively  low  signal  to 
noise  in  these  spectral  regions.  This  may  also  con¬ 
tribute  to  the  disagreement  between  the  model  and 
data  in  the  high  transmittance  regions.  Figure  4 
shows  a  comparison  of  the  data  and  model  over  a 
1 -cm-1 -wide  spectral  region.  Note  again  the  excel¬ 
lent  agreement  near  the  line  centers. 


Fig.  4.  Comparison  of  AEDC  hot  cell  data  with 
model  results  over  a  narrow  spectral 
region. 

Effect  of  Doppler  Shift 

Directed  velocity  (as  opposed  to  random,  ther¬ 
mal  velocity)  causes  the  line  center  locations  to 
shift  due  to  the  Doppler  effect.  A  sample  two-zone 
calculation  was  done  to  demonstrate  this.  The  con¬ 
ditions  in  each  zone  are:  1000  K  temperature, 
0.013-atm  total  pressure,  0.001 3-atm  C02  partial 
pressure,  and  a  5.3-cm  path  length.  Thermody¬ 
namic  equilibrium  is  assumed.  The  gas  in  one 
zone  is  moving  at  1 .5  km/sec  towards  the 
observer,  while  the  gas  in  the  other  zone  is  moving 
at  1 .5  km/sec  away  from  the  observer.  This  can  be 
regarded  as  a  simplistic  model  for  the  radial  veloc¬ 
ity  distribution  in  a  rocket  nozzle.  This  velocity  will 
cause  the  line  centers  in  the  zone  moving  towards 
the  observer  to  be  shifted  approximately  0.01  cm-1 
to  the  blue,  while  the  line  centers  in  the  zone  mov¬ 


ing  away  from  the  observer  will  be  shifted  a  like 
amount  to  the  red.  Note  that  the  total  gas  path  is 
the  same  as  that  used  in  the  Phillips/AEDC  experi¬ 
ment  above.  This  allows  the  unshifted  and  shifted 
computed  spectra  to  be  compared.  These  compar¬ 
isons  appear  in  Figs.  5  and  6.  Figure  5  compares 
the  computed  transmittance  with  and  without  the 
Doppler  shift  over  the  same  1 -cm-1 -wide  spectral 
region  as  seen  in  Fig.  4.  Figure  6  shows  the  com¬ 
parison  of  the  same  results  over  a  narrower  spec¬ 
tral  region.  The  splitting  of  the  spectral  lines  due  to 
the  Doppler  shift  is  readily  apparent. 


Fig.  5.  Comparison  of  model  results  with  and  with¬ 
out  Doppler  shift. 


Fig.  6.  Comparison  of  model  results  with  and  with¬ 
out  Doppler  shift  over  a  narrow  spectral 
region. 


Vibrational  Excitation  Behind  a  Normal  Shock 

A  simplistic  treatment  of  the  vibrational  excita¬ 
tion  behind  a  normal  shock  is  seen  in  Fig.  7.  The 
conditions  for  these  calculations  are  shown  in 
Table  1 .  The  vibrational  populations  are  assumed 
to  be  initially  frozen  at  the  upstream  condition  when 
they  first  cross  the  shock.  The  vibrational  popula¬ 
tions  then  approach  equilibrium  with  the  rotational 
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states.  No  explicit  relaxation  calculation  was  per¬ 
formed.  The  vibrational  temperatures  in  Table  1  are 
assumed  to  exist  at  some  location  downstream  of 
the  shock.  The  rotational  temperature  is  assumed 
to  be  in  equilibrium  with  the  translational  tempera¬ 
ture  at  all  conditions.  The  ambient  gas  is  assumed 
to  be  air,  containing  30  ppm  C02.  The  path  length 
is  held  fixed  at  8  cm.  The  three  vibrational  modes  of 
C02  (bending,  symmetric  stretch,  and  asymmetric 
stretch)  were  assumed  to  be  in  equilibrium.  Future 
work  will  use  both  multitemperature  and  state-to- 
state  approaches  for  this  flow  field. 


Fig.  7.  Vibrational  relaxation  behind  a  normal  shock. 


Table  1.  Conditions  for  Normal  Shock  Calculations 


Parameter 

Upstream 

Conditions 

Downstream  Conditions 

Station  Number 

0 

1 

2 

3 

4 

Mach  Number 

9.8 

0.29 

0.29 

0.29 

0.29 

Pressure,  atm 

0.01 

1.14 

1.14 

1.14 

1.14 

Tot,  K 

300 

3000 

3000 

3000 

3000 

TVjb>  K 

300 

300 

1000 

2000 

3000 

The  results  of  these  calculations  are  shown  in 
Fig.  7.  The  results  have  been  normalized  so  that  the 
change  in  vibrational  structure  can  be  seen  in  a  sin¬ 
gle  plot.  The  increasing  vibrational  excitation  as  the 
vibrational  temperature  approaches  the  rotational 
temperature  is  readily  apparent.  The  drastic  change 
in  the  band  profile  as  equilibrium  is  approached 
should  be  experimentally  observable.  More  detailed 
flow-field  and  relaxation  calculations  would  be 
required  to  compute  physically  realistic  spectra. 

Summary  and  Recommendations 

A  previously  developed  line-by-line  spectral 
model  has  been  extended  to  use  the  HITRAN  and 


HITEMP  line  strength  data  bases.  The  new  spec¬ 
tral  model  computes  the  line-by-line  spectral  radi¬ 
ance  and  transmittance  of  homogeneous  or  inho¬ 
mogeneous  gas  paths,  including  the  effects  of 
Doppler  shifts  due  to  directed  gas  velocity  and 
vibrational  nonequilibrium.  The  nonequilibrium 
vibrational  state  populations  can  either  be  speci¬ 
fied  directly  (perhaps  by  a  state-specific  chemistry 
model,  or  by  imposing  rotational  and  vibrational 
temperatures  for  each  molecular  constituent). 
Model  results  were  compared  to  benchmark  labo¬ 
ratory  data  for  hot  and  hot-through-cold  gas  paths 
of  C02  at  high  and  low  spectral  resolution.  The 
model  compared  favorably  to  the  data  in  most 
cases.  However,  it  does  appear  that  the  treatment 
of  the  line  tails  is  inadequate  to  capture  some  fea¬ 
tures.  This  deficiency  will  be  addressed  in  later  ver¬ 
sions  of  the  model.  The  effect  of  Doppler  shifts  was 
demonstrated  in  a  sample  calculation  with  directed 
velocity  components  on  the  order  of  that  expected 
in  rocket  nozzle  exhausts.  Line  splitting  was  clearly 
visible  in  the  results.  A  simplistic  model  of  vibra¬ 
tional  excitation  behind  a  normal  shock  was  pre¬ 
sented,  demonstrating  the  ability  of  the  model  to 
use  a  multitemperature  description  of  the  vibra¬ 
tional/rotational  populations.  The  ability  of  the 
spectral  model  to  interface  with  state-to-state 
vibrational  chemistry  models  will  be  presented  in  a 
later  paper. 
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